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First application of the Trojan Horse Method with a Radioactive Ion Beam: study of 
the 18 F(p, o) ls O reaction at astrophysical energies. 
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Measurement of nuclear cross sections at astrophysical energies involving unstable species is one 
of the most challenging tasks in experimental nuclear physics. The use of indirect methods is often 
unavoidable in this scenario. In this paper the Trojan Horse Method is applied for the first time to 
a radioactive ion beam induced reaction studying the 18 F(p,a) 15 0 process at low energies relevant 
to astrophysics via the three body reaction 2 H( 18 F, a ls O)n. The knowledge of the ls F(p, a) ls O 
reaction rate is crucial to understand the nova explosion phenomena. The cross section of this 
reaction is characterized by the presence of several resonances in 19 Ne and possibly interference 
effects among them. The results reported in Literature are not satisfactory and new investigations of 
the 18 F(p,q) 15 0 reaction cross section will be useful. In the present work the spin-parity assignments 
of relevant levels have been discussed and the astrophysical S'-factor has been extracted considering 
also interference effects. 

PACS numbers: 26.50+x, 25.70.Hi, 29.38.-c 


In many astrophysical scenarios radioactive nuclei are 
fundamental in two aspects: directly because of their role 
in the nucleosynthesis path and indirectly as information 
carriers that allow to develop and check astrophysical 
models. 

Unfortunately, the difficulties typical of measurements 
of nuclear reaction cross sections at astrophysical ener¬ 
gies greatly increase when using Radioactive Ion Beams 
(RIBs). This is mainly due to the available intensities 
of these beams, at least three orders of magnitude lower 
than the stable ones. In this framework, the development 
and the use of reliable indirect methods to measure these 
cross sections becomes even more important than in the 
case of experiment with stable beams. 

The Trojan Horse Method (THM) has been widely de¬ 
veloped and exploited over the past two decades to study 
bare nuclear cross sections at very low energy and it is 
now routinely used to study reactions between charged 
particles even at zero energy. A detailed description of 
the method can be found elsewhere [Ml and goes be- 
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yond the scope of this paper: a review of the method 
can be found in fl3| . Here we present a pioneering work 
where the THM has been applied for the first time to a 
reaction induced by a RIB, namely the 18 F(p,a) 15 0 pro¬ 
cess at nova energies. 

The 7 -ray emission following the nova explosion is dom¬ 
inated by the 511 keV energy line, coming from the an¬ 
nihilation of positrons produced by the decay of radioac¬ 
tive nuclei. Among them, 1S F is especially important be¬ 
cause of its expected abundance in the Nova environment 
and because of its lifetime, that matches well with the 
timescale for the Nova ejecta to become transparent to 
7 -ray emission. In order 7 -ray astronomy to be helpful in 
understanding the Nova explosion phenomena, it is then 
crucial to know the rate of the nuclear reactions produc¬ 
ing and destroying 18 F. Indeed, at relevant temperatures 
(Tg ss 0.2 — 0.4), which corresponds to Gamow windows 
energies in the center-of-mass of K c m . ss 100 — 400 keV, 
the 18 F(p,a) 15 0 reaction is expected to dominate over 
the 18 F(p, 7 ) by roughly a factor of 1000 and it is the 
uncertainty in this reaction rate that gives the main nu¬ 
clear contribution to the overall uncertainty in the final 
abundance of 18 F jl4|, UB|- In the energy range of inter¬ 
est the cross section of this reaction is dominated by the 
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FIG. 1: Q-value spectra. The Q-value spectrum for the 
2 H( 18 F, a 15 0)n reaction (solid line) is overlaid to the total 
Q-value spectrum (dashed line). The arrow represents the 
theoretical Q-value (0.658 MeV) for the reaction at hand. 
The line is a Gaussian fit with y= 0.668 MeV and <r=0.322 
MeV, in agreement with the expected theoretical value and 
the experimental resolution. 


contribution from states in the 19 Ne compound nucleus 
around the 18 F+p threshold. Especially important is the 
knowledge of the width and spin-parity of these resonance 
states, in order to verify possible interference effects af¬ 
fecting the rate. However, in spite of many and long 
lasting attempts to infer information on the 18 F(p,a) 15 0 
reaction by both direct and indirect measurements (see 
e.g. }16h32| ). the situation is still not satisfactory. In 
order to get new complementary information on this pro¬ 
cess, the 18 F(p,a) 15 0 reaction has been studied in inverse 
kinematics by applying the THM to the three body re¬ 
action 2 H( 18 F, a 15 0)n. If one succeeds in selecting the 
phase space region where this reaction proceeds through 
a quasifree reaction mechanism, then its cross section can 
be factorized [Ml- The bare nucleus cross section of the 
18 F(p,a) 15 0 reaction can hence be easily deduced even 
at very low energies overcoming the suppression effects 
coming from coulombian and/or centrifugal barriers pen¬ 
etrability factors and the uncertainties due to the poor 
knowledge of spectroscopic information. 

The fundamental step in applying this method is the se¬ 
lection of the events that proceed through quasifree re¬ 
action mechanism among those coming from all other 
possible mechanisms. To this aim, two out of the three 
reaction products must be detected with typical energy 
and angular resolution of 1% and 0.1°-0.2° respectively. 
In particular, some reconstructed kinematical variables 
(like relative energies of the outgoing particles) are ex- 
pecially sensitive to the angular resolution of the exper¬ 
imental setup. Given the characteristics of the radioac¬ 
tive beam used in the present experiment and to achieve 
the angular resolution required by the THM, the exper¬ 
imental setup 0 was designed so that it was possible 
to track the beam particles event by event and to detect 
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FIG. 2: Momentum distribution for the p-n intercluster mo¬ 
tion in deuteron. The solid line is the Hulthen function in 
momentum space. 

ejectiles going in a large solid angle. The beam tracking 
system, based on a pair of Parallel Plate Avalanche Coun¬ 
ters (PPACs), allowed for an angular resolution of 0.14°. 
The two ejectiles were detected by using two planes of 
position sensitive detectors covering the angular ranges 
2°-ll° and 11°-31° with respect to the geometrical cen¬ 
ter of the target. In particular a pair of double-sided 
multistrip silicon detectors (DSSSD), was used to detect 
the heaviest ejecta. The lighter ones were detected by 
using the ASTRHO (A Silicon Array for TRojan HOrse) 
modular system of INFN Laboratori Nazionali del Sud 
equipped with 8 bidimensional position sensitive detec¬ 
tors (BPSD, 45x45 mm 2 , 500 ^m thick, spatial resolu¬ 
tion 1 mm) similar to those described in [33l . (t 3 |. For 
each detector the x and y position and the energy of the 
impinging particles were recorded together with the time 
of flight. A standard energy calibration procedure was 
used for all detectors while the x-y position coordinates 
given by the BPSD were calibrated by placing grids in 
front of them. 

The emission angles of the ejectile tracks with respect 
to the measured beam track were reconstructed with an 
overall angular resolution of 0.2°. 

The 1S F beam was produced via the 18 0(p,n) 18 F reac¬ 
tion using the CNS Radioactive Ion Beam (CRIB) sepa¬ 
rator of the Center for Nuclear Study (CNS) of the Uni¬ 
versity of Tokyo, installed at RIKEN accelerator facility 
in Wako, Japan [35j,[3(|. The characteristics of the 18 F 
beam during the experiment were the following: beam 
energy peaked at 47.9 MeV (FWHM 1.9 MeV), maxi¬ 
mum beam intensity of 2xl0 6 pps and purity better than 
98%. The intensity was higher than 5x10 s pps through¬ 
out the measurement. The secondary target was made of 
a thin (typically 150 /ig/cm 2 ) CD 2 foil. The chosen beam 
energy fulfills the prescriptions for the validity of the Im¬ 
pulse Approximation |37( : (a) the wavelength associated 
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to the entry channel is smaller than the nuclear radius of 
deuteron (A = 1.54 fm), in order to maximize the proba¬ 
bility that quasifree reaction mechanism occurs; (b) the 
incident center-of-mass energy (4.82 MeV) is higher than 
the binding energy of deuteron. 

The used setup assures that the coincidence efficiency 
is independent from the excitation energy. A detailed 
Monte Carlo simulation was performed for the used 
geometry and kinematics assuming that the a+ 15 0 
breakup is isotropic in the center-of-mass system. The 
actual beam energy profile, the reconstructed beam 
tracks distribution and the measured beam position on 
the target were used in the simulation. The result of the 
simulation was used to calibrate the solid angle cover¬ 
age of the detectors. In order to select events coming 
from the reaction of interest, several cuts were applied 
on the data set. In particular it was requested that only 
one of BPSD and one strip of one DSSSDs fired with 
multiplicity equal to 1 on each detector and that the 
time of flight of the two particle was within the corre¬ 
lation time window. The events coming from reactions 
of 18 F on carbon were easily removed thanks to the dif¬ 
ference in phase space distribution of these events with 
respect to those induced on deuterium. Moreover, study¬ 
ing the bidimensional spectra of the relative energy of 
one pair out of three outgoing particles versus that of 
other pairs (e.g. E a _iso vs E a _„), it was possible to 
separte the phase space region spanned by the events 
coming from the 2 H( 18 F, a 15 0)n reaction that partially 
overlaps with those coming from other possible reactions 
on deuterium, namely 2 H( 18 F, a 15 N)p, 2 H( 18 F,p 18 0)p 
and 2 H( 18 F,p 18 F)n. Fig. |T] shows the Q-value spec¬ 
trum that was obtained by applying this event selection 
procedure (solid histogram) overlaid to the total Q-value 
spectrum (dashed histogram). The line shown in this fig¬ 
ure is a Gaussian fit with p= 0.668 MeV, to be compared 
with the theoretical Q-value of 0.658 MeV, and er=0.322 
MeV, in agreement with the overall experimental reso¬ 
lution. The good agreement between the experimental 
and the theoretical value gives confidence in the correct 
identification of the events coming from the reaction of 
interest. 

In order to apply the THM, it is essential to select the 
events coming from the quasifree reaction mechanism in 
the 2 H( 18 F, a 15 0)n reaction among other reaction mech¬ 
anism. The strongest evidence of the predominance of 
the quasifree mechanism is given by the shape of the 
momentum distribution for the p-n intercluster motion 
in deuteron [§.]. Data are shown in Fig. [2] by black dots. 
The solid line in Fig. [2]represents the Hulthen function in 
momentum space with standard parameters values [38t | . 
The fair agreement between data and theoretical curve 
give confidence in the selection of events coming from 
the quasifree reaction channel and allows to use a Plane 
Wave Impulse Approximation in the data analysis [38| . 
Indeed, it is well known that distortion effects influence 
the behaviour of the momentum distribution at higher 
values of momentum. Consequently, the rest of the anal- 



FIG. 3: The nuclear cross section spectrum in function of 
the p- 18 F cm energy for the events that pass the conditions 
described in the text. The blue vertical line shows the position 
of the threshold for the ls F+p reaction (Eth = 6.41 MeV). 
The red dashed lines represent the gaussians used for fitting 
the data as explained in the text. The numbers above the 
arrows represents the peak positions in 19 Ne excitation energy 
obtained from the fitting procedure. 

ysis was done choosing the events having a spectator mo¬ 
mentum lower than 60 MeV/c. Assuming that the events 
selected according to the previous procedure come from 
the quasifree contribution to the reaction yield, the cross 
section of the 18 F(p,a) 15 0 process can be deduced divid¬ 
ing the yield of the 2 H( 18 F, a 15 0)n reaction by a kine- 
matical factor and the momentum distribution of p and n 
inside the deuteron [Ml- The nuclear cross section for 
18 F(p,a) 15 0 reaction, unaffected bv suppression effects 
due to Coulomb and centrifugal [12J barriers, measured 
by THM down to zero (and even negative) p+ 18 F rel¬ 
ative energy is shown in Fig. [3] The obtained spectra 
were analyzed using a least-square fit of multiple Gaus¬ 
sian. The sigma of the Gaussians was fixed to 53 keV 
based on the energy resolution calculated for the present 
experiment. The excitation energies coming from the fit 
together with the error on the peak position are listed 
in Table I. For comparison energies and J’ 1 ’ coming from 
Ref. [ 31 ], |H, HI, HI are also reported in Table I. 

Further analysis have been done by considering only 
the energy region of interest for astrophysical purposes, 
i.e. the observed excited states at energies 6255, 6460, 
6537, 6755, 6968 and 7075 keV. 

Though the cross section of the 18 F(p,a) 15 0 reaction 
was indirectly measured by THM in the center of 
mass angular range 70° < 6 cm <120°, the statistics of 
the present data do not allow to extract the angular 
distributions of the populated resonances and hence 
a self-assignment of the values is not possible in 
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E c .m. (keV) 

E (keV) 

Ref. [39] Ref. [40] 

Ref. [311 

Ref. [32] 

(present work) 

(present work) 

C.D. Nesaraja et al. D.R. Tilley et al. 

A.S. Adekola et al. 

A.M. Laid et al. 

-574±17 

5837 

5832±9 





6013±7 (3/2,1/2") 


6014±2 (3/2") 

-341±16 

6070 


6089±2 

6072±2 (3/2+,5/2”) 



6092±8 


6097±3 (7/2,9/2)+ 



6149±20 


6132±3 (3/2+ ,5/2“) 

-156±18 

6255 

6288±7 

6289±2 (l/2 + ,3/2 + ,l/2“) 

6289±3 (5/2,11/2)" 


6419±6 (3/2+) 6419±6 (3/2") 6416±3 (3/2",5/2+) 

(6422)±30 (11/2+) 

6437±9 (l/2“) 6437±9 6440±3 (11/2+) 

6449±7 (3/2+) 

49±14 6460 (3/2+, 5/2") 6459±3 (5/2") 


126±15 

6537 (7/2+,9/2+) 

(6504)±30 (7/2+) 

(6542)±30 (9/2+) 






6698±6 (5/2+) 



6700±3 

344±18 

6755 (3/2") 

6741±6 (3/2") 

6742±7 (3/2,1/2-) 

6747±5 (3/2“) 

6742±2 (3/2") 



(6841)±30 (3/2“) 
6861±6 (7/2+) 

6861±7 


6862±2 (7/2") 

556±19 

6968 (5/2+) 

(6939)±30 (l/2“) 

(7054) ±30 (5/2+) 




664±10 

7075 (3/2+) 

7075.7±16 (3/2+) 

7067±9 

7089±5 (3/2+) 



924±11 

7335 


7326±15 



1056±13 

7467 

7420±14 (7/2+) 
7500±9 (5/2+) 

7531±15 

7431 



TABLE I: Summary of the 19 Ne resonance parameters in the energy range explored by the experiment compared to results 
from other works. The energies of Ref. [39j] reported in parenthesis have been not measured. 


this case. However, as the THM data are not affected 
by suppression effects coming from the Coulomb and 
centrifugal barrier 01, the population of the excited 
states in the compound nucleus is linked to the J 77 of the 
levels (4lj| . Hence, for each populated 19 Ne excited states 
observed in the present experiment a certain J 7r value 
has been assumed by comparison with those available 
in Literature, as discussed below, and it is reported in 
parenthesis in Table I. Using these assumptions, data 
from each resonance have been integrated over the full 
angular range by means the corresponding Legendre 
polynomial to obtain the contribution to the THM 
cross section. Finally, each contribution to the THM 
cross section has been multiplied for the corresponding 


penetrability factor of the centrifugal barrier to get the 
astrophysical S'-factor of the reaction. As THM cannot 
provide for the absolute normalization of results, the 
obtained spectrum for the S'-factor was then normalized 
to the results of direct measurements of the well known 
resonance at 7075 keV. 

From Literature, the J 7r values of the two resonances at 
energies 7075 keV and 6755 keV are well known to be 
3/2 + and 3/2“, respectively, and these values have been 
assumed throughout the present analysis. 

The resonance at E = 6968 keV has not been observed 
in previous measurements. Nonetheless, following 
Ref. [39], two nearby levels are predicted at 6939 keV 
(J 77 = 1/2-) and at 7054 keV (J 77 = 5/2+). The 
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resonance at E = 6968 keV could then be associated 
with these predicted levels. In this case, fixing a value 
of J 77 = 1/2 - results in a cross section value that is 
unlikely bigger than that of the resonance at E = 7075 
keV. Moreover, the fact that the resonance at E = 6968 
keV is not evident in direct measurements suggests 
instead a relatively high l assignment. So J 77 = 5/2+ 
(l = 2) has been assumed for this level in the present 
work for the calculation of the astrophysical S-factor, 
though we cannot exclude higher assignments of the 
value of l for this resonance. 

With a similar reasoning, the 19 Ne state at E = 6537 
keV has been associated with the calculated levels at 
6504 keV (J 77 = 7/2+) or 6542 keV (J 77 = 9/2+) both 
with l = 2 [39j]. Both values of the resonance parameters 
result in a strong suppression of the cross section due to 
the centrifugal barrier penetrability and hence do not 
influence the final value of the astrophysical S'-factor. 
Just for computational purposes the spin-parity assign¬ 
ment for this level was fixed to be 7/2+. 

The situation regarding the lower energy region is less 
clear. In the present experiment only a level at 6460 
keV is observed. Above the proton threshold several 
resonances are foreseen by comparing 19 Ne to its mirror 
nucleus 19 F [31] and some of them have been observed 
OlIlEl- In particular, a resonance at E = 6459 keV 
has been observed in [33 and it has been interpreted 
as part of a triplet of states with possible spins and 
parities (3/2“ or 5/2+) at 6416 keV, (11/2+) at 6440 
keV and (5/2 _ ) at 6459 keV. Calculations presented in 
[39| attribute to four states in the same energy region 
the spin-parity values of 1/2“ (6419 keV), 3/2+ (6422 
keV and 6449 keV) and 11/2+ at the unobserved state 
at 6422 keV. Though only the level at 6449 keV and 
that at 6459 keV are within the fit error for the 6460 
keV peak observed in the present work, calculations 
of the contribution to the total cross section due to 
this very level have been performed assuming all of the 
J 77 values mentioned above. If the spin-parity value is 
fixed to be 11/2+ or 5/2+ this contribution is strongly 
suppressed by the centrifugal barrier penetrability factor 
and hence these J 77 assignments are rejected. On the 
other hand, there is no reason to rule out the other 
values of J 77 for this level, namely 1/2“, 3/2“, 5/2“ and 
3/2+. Calculations showed that the differences on the 
contribution of the 6460 keV level to the astrophysical 
S'-factor for spin-parity assignment 1/2“, 3/2“ and 
5/2“ are negligeable within the errors. To conclude the 
discussion on the 6460 keV level, it is worth noting that 
possible interference effects in THM are not calculated 
but are already contained in the data [H, [4l|. So the 
other possible assignment 3/2+ to the 6460 keV level 
will automatically take into account interference effects, 
if any. 

Finally, in the subthreshold region, the excited state ob¬ 
served here at E = 6255 keV was assigned a J 77 = 11/2 - , 
as already proposed in Ref. j.32[ : any other assignment 
on the base of the nearby levels observed in |3j], (32] 


results in an unrealistically huge contribution of this 
level to the total astrophysical S'-factor. 

In Fig. [4] the results obtained in this work for the astro- 
physical S-factor are presented: assuming J 77 = 3/2+ 
for the 6460 keV state the result is reported as full dots 
while the J 77 = 5/2“ assumption for the same level gives 
the astrophysical S-factor shown as open dots. In Fig. 
U the experimental points are also compared with the 
calculations for the as trop hysical S-factor presented in 
Fig. 3 and 4 of Ref. [3C| smeared at the experimental 
resolution obtained in this work. In particular, the 
solid lines represent the upper and lower limits for an 
R-matrix calculation where the interference among the 
three states with J 77 = 3/2+ at energies 6419, 6449, 7075 
keV has been considered. The dashed lines represent the 
same limits for the case with interference between the 
two states at E = 6449 keV and E = 7075 keV only. In 
this latter calculation the authors of Ref. [30| attributed 
the value J 77 =3/2“ to the E = 6419 keV state in 19 Ne. 
Data from the present experiment have been normalized 
to these calculations imposing that the integral of the 
resonance at 7075 keV is the same in direct and THM 
data. 

In the energy region below 100 keV, depending on the 
J 77 assignment chosen for the 6460 keV level, the present 
data either agree fairly well with the region given by the 
dashed lines in Fig. []] (J 77 = 3/2+, full dots) or become 
much lower (J 77 = 5/2“, open dots) than any previous 
result. In both cases the data obtained in this work 
seems to exclude the existence of three interfering states 
having J 77 = 3/2+ represented by the solid lines Fig. 0] 
In conclusion, the THM was applied for the first time 
to study a reaction induced by a radioactive ion beam. 
Even with the use of indirect methods like THM, the 
measurement of cross sections of interest for nuclear 
astrophysics remains one of the most difficult tasks in 
nuclear physics, because the low radioactive ion beam 
intensity adds on the top of the low cross sections typical 
of astrophysical nuclear processes. 

The THM data have been used to obtain the nuclear 
cross section for the 18 F(p,a) 15 0 reaction and, by com¬ 
parison with pieces of information present in Literature, 
to infer information about the J 77 of the 19 Ne nucleus 
excited states. From this it was possible to extract the 
astrophysical S-factor for the 18 F(p,a) 15 0 process. In 
particular resonances in 19 Ne at energies 6255, 6460, 
6537, 6755, 6968, 7075 keV have been observed and 
studied, as they mostly influence the energy region of 
interest for the nova phenomena. A value of J 77 = 5/2+ 
has been assigned to the excited state at E = 6968 keV. 
For the E = 6537 keV both J 77 = 7/2+ or 9/2+ are 
compatible with the present data as this contribution 
to the astrophysical S'-factor is negligeable. In the 
subthreshold region, the excited state at E = 6255 keV 
was assigned a J 77 = 11/2 - , following [32j . 

Finally, a single resonance has been observed just above 
the proton threshold at E = 6460 keV. Different spin 
assignments have been considered, namely J 77 = 3/2+ 
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E (keV) 

FIG. 4: The 18 F(p, a) 15 0 astrophysical S'-factor from the 
present experiment. The full dots are THM experimental data 
with the assumption of J n = 3/2 + for the resonance at E = 
6460 keV, the open ones corresponds to the assumption of 
J 77 = 5/2“ (the difference from this last assumption to the 
other possible value 1/2“ and 3/2“ being negligeable within 
the errors). The solid and dashed lines shown in figure are 
calculations presented and discussed in Ref. [§0j] smeared to 
the present experimental resolution: see text for details. 


(Fig. [4] full dots) and J 77 = l/2“, 3/2“ or 5/2“ (Fig. [4] 
open dots). The comparison between the experimental 
data and the calculations of Ref. [30] seems to exclude 
the presence of two excited states near the proton 
threshold having J 77 = 3/2 + . It is worth noting that, 
if the value J 77 = 3/2 + is assigned to the E = 6460 
keV state, the effects coming from its interference with 
the E = 7075 keV state are implicitly present in the 
THM data. Therefore, once measured unambigously 
the J 77 of this excited state, a straightforward indication 
of the value of the astrophysical S-factor -including 
interference effects- in the energy region of interest for 
Novae could immediately come from THM data. 
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